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Abstract 

The effect of quenched impurities on systems which undergo first-order phase 
transitions is studied within the framework of the g-state Potts model. For 
large q a mapping to the random field Ising model is introduced which pro- 
vides a simple physical explanation of the absence of any latent heat in 2D, 
and suggests that in higher dimensions such systems should exhibit a tricrit- 
ical point with a correlation length exponent related to the exponents of the 
random field model by v = z^rf/(2 — orf — /3rf)- I n 2D we analyze the model 
using finite-size scaling and conformal invariance, and find a continuous tran- 
sition with a magnetic exponent f3/v which varies continuously with q, and a 
weakly varying correlation length exponent v w 1. We find strong evidence 
for the multiscaling of the correlation functions as expected for such random 
systems. 
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Although the effect on the critical behavior of adding quenched bond randomness to 
classical systems whose pure version undergoes a continuous phase transition is well under- 
stood in terms of the Harris criterion |1| , the analogous situation when the pure transition is 
first order is less well studied. Following earlier work of Imry and Wortis ||, Aizenman and 
Wehr [[J and Hui and Berker j|] argued that in 2D any amount of randomness should lead 
to a vanishing of the latent heat. The arguments leading to this conclusion are analogous to 
those used by Imry and Ma || for the absence of any spontaneous magnetization in the ran- 
dom field Ising model (RFIM) for d = 2: the bond randomness couples to the local energy 
density, which is different for the coexisting phases of the pure model, in the same manner 
that the random field couples to the local magnetization of the RFIM. The vanishing of the 
latent heat should be accompanied by a divergent correlation length, and, if so, the question 
arises as to which universality class(es) the corresponding continuous transition belongs. A 
suitable model in which to study this is the g-state Potts model, whose pure version in 2D 
undergoes a first-order transition for q > 4, otherwise being continuous. Chen, Ferrenberg 
and Landau || undertook an extensive Monte Carlo investigation of the case q = 8. In 
addition to confirming the continuous nature of the transition, these authors extracted nu- 
merical values of the critical exponents which appear to be consistent with those of the pure 
2D Ising model. Similar values have also been claimed for the case q = 4, when the pure 
transition is continuous @. This disagrees with the predictions of Ludwig and Cardy |§, 
Ludwig || and Dotsenko et al. who find a new random fixed point for q > 2, based on 



an expansion in powers of q — 2. 

In this Letter we describe the results of an investigation of this model along two different 
lines. First, we show that at large q the model (or rather a model for the interface between the 
disordered and ordered phases) may be mapped onto a corresponding problem for the RFIM. 
This mapping gives a direct and intuitive picture for the arguments of Refs. . Moreover, 
since the RFIM interface is well understood from the renormalization group (RG) perspective 



n|, one may immediately transcribe these results to give quantitative predictions for the 



random bond problem, both in two and in higher dimensions. The second line uses the 



powerful techniques of conformal invariance, combined with finite-size scaling, which should 
apply to such a system with a divergent correlation length. By generalizing the transfer 
matrix formalism of Blote and Nightingale [12] we are able to generate results for arbitrary 
q and reasonably large strip widths. At the same time we have to take account of the 
fact that in random systems the transfer matrices do not commute, and one must discuss 
the spectrum of Lyapunov exponents rather than the eigenvalues. The relation of these to 
critical exponents is indirect, but we have developed a method of making this connection 
based on a cumulant expansion, in which we see explicitly the multiscaling properties of 
the correlation functions discussed theoretically by Ludwig |J. The above method works 



poorly for the thermal exponent. Instead, by finite-size scaling arguments [J3| we are able 
to determine v directly. 

As will become clear, many of our results generalize, but let us for definiteness consider 
a Potts model on the square lattice with degrees of freedom Si taking q values, and a 
reduced hamiltonian —^ijKij5 SiSp where the sum is over nearest neighbor pairs. The 
ferromagnetic couplings Kij are quenched random variables, taking the values K\ and K 2 , 
each with probability |. When (e Kl — l)(e K2 — 1) = q this model is, on average, self- 
dual, and, if the transition is unique, is therefore at its critical point ]TJ|. It is useful to 



parametrize e Kii — 1 = Uij = q2 +Wi \ where Wy = ±w, and w > measures the strength 
of the randomness. The partition function of this model may be mapped onto that of the 



random cluster model []15]| , in which each bond of the lattice is either occupied, when it is 
counted with weight u^, or empty, in which case it is counted with weight 1. The partition 
sum is over all such configurations, in which each connected cluster of sites is weighted by a 
factor q. Let us first consider the pure model, with w = 0. In the limit q — > oo, the sum over 
configurations is dominated by only two: the empty lattice, in which no bonds are occupied, 
which contributes a factor q , where N is the total number of sites, and the full lattice, 
with a weight {^fq) 2N ' , since the number of bonds is 0(2N). All other configurations are 
down by powers of q^ . At the self-dual point, there are therefore two co-existing states with 
identical bulk free energy and different internal energy densities, indicating, as expected, 



that the transition has a non- vanishing reduced latent heat per bond ~ \ In q. For the pure 
model, this analysis may be extended to take into account higher order corrections in q~^, 
with no essential change in the physical picture. Now consider an interface between these 
two phases. For large q, the lowest energy interface is parallel to a lattice direction, say 
the x-axis, and is such that all the bonds with y < some integer are occupied, and those 
above this are empty (or vice versa.) There will also be entropic fluctuations y = h(x) of 
this interface, described by the usual solid-on-solid interfacial hamiltonian, proportional to 
the length of the interface. The interfacial tension for large q is a ~ \ In q, independent of 
the local shape of the interface. This is to be compared with a ~ 2J between the ordered 
phases of a low temperature Ising model with reduced exchange coupling J. 

Now consider the effect of adding bond randomness to the random cluster model. Each 
configuration of the interface will be weighted by an energy J2 X H y <h(x) w(x, y) lng, where 
(x, y) labels bond positions. This may be rewritten, up to a term independent of h(x), as 
\ J2 X ( J2 y <h(x) ~ J2 y >h(x) )w(x, y) lng. Compared with the energy of an RFIM interface with 
spins s(x,y) = ±1 coupled to a reduced random field h(x,y) = ±/i RF , the interfacial models 
are identical with the correspondence J <-> |lng, and h RF <-> |uTng. In addition, the 
imposition of a uniform reduced magnetic field h on the RFIM, which distinguishes between 
the two co-existing phases, is seen to be equivalent to a deviation t = (T — T c )/T c in the 
temperature variable away from the critical self-dual point. Since this couples to the energy 
density we find the correspondence h <-> \t In q. 

Of course, this is strictly valid only as q — > oo. At finite q the g-dependence of cluster 
configurations with more complicated topologies is not simply accounted for by the interfacial 
tension. For the same reason, the mapping is not between bulk configurations of the two 
models. However, it will be argued that certain universal properties are controlled by an RG 
fixed point at infinite q, and for these the mapping should be asymptotically exact. Although 
this has been described in terms of a 2D self-dual model, it should be clear that it is more 
general: lack of self-duality corresponds to a skewness in the distribution of the random 
fields h(x,y), which may be compensated by adding a suitable uniform field (corresponding 
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to a shift in the T c of the Potts model), and, similarly, higher dimensions may be taken into 
account by appropriately replacing ^/q by q 1 ^. 



The RG properties of the interface in the RFIM near d = 2 have been well studied fllT 
When translated into the variables of the random bond Potts model, the flow equations have 
the form 

dw/dl = -(d/2-l)w + Aw 3 + --- (1) 
dilnqY^dl = -(lnq)-\(d -I)- Aw 2 + ■■■) (2) 
dt/dl = t(l + Aw 2 + •••), (3) 

where A > is a non- universal constant. Corrections to these equations are supposed to 
be higher order in w and in q" 1 ^ 2 . The reader who is concerned that the number of Potts 
states q is allowed to flow in these equations is reminded that q is merely a parameter of 
the random cluster model. The RG flows for d > 2 and the consequent phase diagram are 
shown in Fig. |1| In the pure models, for q > some qi{d) (low T in the RFIM), there is phase 
coexistence with a non- vanishing latent heat (spontaneous magnetization), controlled by a 
fixed point at infinite q (T — 0). For d > 2 this persists into the shaded region, bounded 
by a line of tri critical points where the latent heat, vanishes. The universal behavior along 
this line is controlled by the fixed point R at w = 0((d — 2) 1 / 2 ) and infinite q. Using the 
correspondence t ■ T^p, we conclude that the thermal eigenvalue v~ l of the random 

bond problem is related to the RG eigenvalues and exponents of the RFIM by the scaling 
relations v~ l = y h — 9 = (2 — a RF — j3 RF ) / z/ RF , where —6 is the eigenvalue which controls the 



irrelevance of Trf arid the consequent violation of hyperscaling |TT| . In the same manner, 
it may be shown that the latent heat vanishes as (w c — w)@ RF as the line Rq^ is approached 
from below. Of course, these relations have been established only close to d — 2, but, if 
the topology of the RG flows does not change, they should hold also in three and higher 
dimensions. 

Above the line Rq2, the flows go to large w beyond the validity of ([j]-|3]). In addition, the 
renormalized interfacial tension flows to zero and the mapping between the models breaks 



down as domains of different topologies proliferate. However, for infinite q the mapping 
remains exact and the flows go to infinite w. This cannot happen for finite q since this is 
the percolation limit KyjK-i = 0, at which w -1 is relevant ||16| . There must therefore exist 
another line of stable fixed points emerging from P\, which control the universal continuous 
transition for large, but finite, values of w and q. It is tempting to conjecture, as indicated 
by the dashed line in Fig. [I], that this connnects on to that found by expansion in powers 
of q — qi ||, where qi is the point where the exponent a of the pure model changes sign 
0]. Our analysis indicates that, at least for d = 2, this is the case. In 2D (when q\ = 2, 
q2 = 4), the shaded region collapses, and for any non-zero w the renormalized interfacial 
tension, and thus the latent heat, vanish. The flows should be towards the line P\qx, with a 
cross-over length which, from Eq. (TJ), has the form £x ~ ( Al 2Aw ' 1 anc i therefore may become 
very large for weak randomness. 

We now turn to our numerical results for d = 2. These will be described in detail 



elsewhere ]I7||. As shown by Blote and Nightingale JT2| the transfer matrix for the pure q- 
state Potts model in a strip of width L may be constructed in a basis in which q enters only 
as a continuous parameter. We have generalized this to the case when the bond strengths 
Uij are quenched random variables, and the transfer matrices % therefore depend on the 
row labels i. The size of the transfer matrices grows like 4 L , independent of q, making this 
a practicable method for larger q. 

Starting with some suitable initial vector v , the leading Lyapunov exponent is given 
by A° = linv^oo — In [ ] (112= i ^) v o||- Higher exponents are found by iterating a set of 
vectors (vj)^ =0 , where a given Vj is orthogonalized to the set (vj)|~Q after each multiplication 
by the transfer matrix. The average free energy per site is then fi = — -^A°. For a system 
exhibiting a first-order transition with a bulk correlation length £ we expect [[12]] ~ 
/«, + 0(L- d e- L/ Z), so that X(L) = \n[f L - /«,] + rflnL ~ const - L/C In Fig. | we 
show this for various values of q and the randomness strength R = KijK\. We see that 
the randomness changes the transition for q = 8 into one with an apparently diverging 
correlation length. In such a case the amplitude of the finite-size correction has the form 
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IK /l ~ foo — vrc'/6L 2 , where d is the effective central charge (since fi is the quenched 
free energy, in a replica formalism this is the derivative of the central charge c(n) wrt the 
number of replicas at n — 0). The value of d was determined by making parabolic least- 
squares of fi versus 1/L 2 PD| . We found the optimum trade-off between statistical errors 
and a reasonable computation time by taking strips of length m = 10 5 and averaging fi 
over 100 independent realizations of the randomness for 1 < L < 8, and 3 realizations for 
9 < L < 12. Data collection was made for each 200 multiplications by %, and the first 2,000 
iterations of each run were discarded in order to eliminate transients. The parabolic fits 
were made by including the data points for Lq < L < 12, where Lq must be chosen large 
enough to justify the finite-size scaling form, and small enough to minimize error bars. From 
the special cases of the Ising model and percolation it appeared that L = 3 is optimal. 

For the random bond Ising model (q = 2) with R = 2 we found d = 0.495 ± 0.006, in 
agreement with de Queiroz' result |2(J d = 0.498 ± 0.003 using the spin basis, and with the 



expected value d = |. For q = 3 and R = 2, our result d = 0.799 ± 0.006 is unable to 
distinguish between the pure value of | and that of d ~ 0.8025 obtained in Ref. || by an 
expansion in q - 2. For q = 4 the results d = 1.003 ± 0.006 for R = 2 and d = 1.010 ± 0.022 
for R = 10 are consistent with each other and the pure value c = 1, but for larger q the R = 2 
results appear to saturate while those for R = 10 show a gradual increase: d = 1.517 ±0.025 
for q = 8 and d = 3.003 ± 0.031 for q = 64. Similar values have recently been reported by 
Picco [^TJ . However, it should be pointed out that these are also very close to those expected 



at the percolation point R —>■ oo. For then the replicated model is the Potts model with 
q n states, so that d = (d/dn)cp otts (q n )\ n=0 = (5v / 3/47r) lng ps 0.689 lng; this is confirmed 
by our transfer matrix calculations fl7| . These are also remarkably close to the pure values 
for 2 < q < 4. We conclude that measurement of the effective central charge does not 
distinguish well between pure, percolative and non-trivial random behavior. 

We therefore turned to the magnetic exponent x\ = (3 jv. Although this may be deter- 
mined in principle by adding a ghost site [0, it proved more efficient to exploit duality 
and relate the spin-spin correlation function G(m) on the strip to the free energy in the 



presence of a seam of frustrated bonds. Details of this relation and how to implement it in 
the transfer matrix will be given elsewhere [T7||. In pure systems, according to conformal 



symmetry |^2| G{m) decays along the strip as exp(— 2nxim/L), so that the difference A/l 
in the free energy per site with and without the seam behaves as 2irx\/L 2 . However, in the 
random system this difference typically has fluctuations which are 0(m -1 / 2 ). This has the 
consequence that, while Afi ~ -^lnG is self-averaging, G is not [^3[. In fact, as shown by 



Ludwig ||, the moments G(m) N (where the overline denotes the quenched average) exhibit 
multiscaling, that is, they scale with dimensions xn which are not, in general, linear in N. 
Since conformal symmetry assumes translational invariance, it refers only to such averaged 
quantities. In practice, we can avoid the lack of self-averaging by performing a cumulant 
expansion 



lnG iV = MnG + iiV i (lnG'-lnG) 2 + ---, (4) 

where the higher cumulants may be directly extracted from the statistical fluctuations in 
Afi. For values of q and R which are not too large, this expansion appears to converge well, 
keeping the first 3 or 4 cumulants. The fact that the higher cumulants are non-zero implies 
multiscaling. Our values for x\ are shown in Fig. |^. For 2 < q < 3 and R = 2 they are in 
perfect agreement with the predictions of the (q — 2)-expansion of Refs. PJlO|. For larger q 



the results appear to increase smoothly with q, with X\(q = 8) = 0.142 ± 0.001. Thereafter 
the cumulant expansion begins to break down. 

Although the thermal exponent v should be related in a similar manner to the first 
gap — A° L in the Lyapunov spectrum, for reasons we do not understand this yields 



results which, if taken literally, appear to violate the bound v > 1 p4|. Instead we have 
measured v directly by finite-size scaling of the magnetic correlation length away from T c , 
using phenomenological RG methods [|13j. For q = 8, as shown in Fig. |], we find clear 
evidence for fixed points at R — 0, with v = | (xt = 0), and a random fixed point with 
i/r = 1.01 ± 0.02. At q = 3 our results are consistent with the perturbative value z/r 1.02, 
and for q = 64 we find z/r, = 1.02 ± 0.03. 
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In summary, we have computed the exponents of the 2D random bond Potts model 
and shown that while the thermal exponent z/r is consistent, within error bars, with both 
the pure Ising value, and with the results of the (q — 2)-expansion, the magnetic exponent 
varies continuously with q in a manner which agrees with this expansion in the region it is 
expected to be valid, and with a value at q = 8 which is quite different from the Ising value 
of ~. This is in sharp disagreement with the Monte Carlo results of Ref. |J. One possible 
reason is that these authors use a non-standard definition of the order parameter which, 
while it scales in the same way as the usual one in the pure case, may not when multiscaling 
is present. Another is that our very long strips are able to accommodate large regions in 
which all q values of the order parameter are realized, and this may not typically be the 
case in the square geometries of Ref. || . We have also given a mapping of the random bond 
problem to the RFIM which suggests that for d > 2 such systems should exhibit a tricritical 
point whose thermal exponent is related to those of the RFIM. This picture is quite generic, 
and, since there are many real 3D systems which undergo first-order transitions, it would 
be interesting to re-examine the effect of random impurities in such cases. 

The authors acknowledge useful discussions with J. Chalker, A. Ludwig and S. de 
Queiroz. This research was supported in part by the Engineering and Physical Sciences 
Research Council under Grant GR/J78327, and by the National Science Foundation under 
Grant PHY94-07194. 
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FIGURES 
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FIG. 1. Schematic phase diagram in the critical surface for d > 2. q increases to the left and 
w is the disorder strength, with -P1-P2 being the percolation limit. RG flows are indicated. The 
latent heat is non-vanishing within the shaded region, and elsewhere the transition is continuous, 
controlled by the line of fixed points P\q\. 
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FIG. 2. Plots of A(L), normalized to A(l) = 1, showing that bond randomness renders the 
phase transition second order (see text). 
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FIG. 3. Magnetic exponent x\ = {3/v as a function of q. 
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FIG. 4. Values of ^ extracted from phenomenological RG for different strip widths and q = 8. 
There is a fixed point at R « 8 with f « 1 (see inset). Error bars are less than the symbol size. 
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